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Abstract 

Modelling of inspection data for large scale physical systems is critical to assessment of their integrity. We 
present a general method for inference about system state and associated model variance structure from 
spatially distributed time series which are typically short, irregular, incomplete and not directly observable. 
Bayes linear analysis simplifies parameter estimation and avoids often-unrealistic distributional assumptions. 
Second-order exchangeability judgements facilitate variance learning for sparse inspection time-series. The 
model is applied to inspection data for minimum wall thickness from corroding pipe-work networks on a 
full-scale offshore platform, and shown to give materially different forecasts of remnant life compared to an 
equivalent model neglecting variance learning. 

Keywords: Bayes linear, exchangeability, variance learning, corrosion, dynamic linear model, Mahalanobis 
distance 



1. Introduction 

Inspection and maintenance assures the integrity of physical systems subject to degradation (e.g. to cor- 
rosion and fouling in time). Inspection and maintenance is typically costly, requiring careful allocation of 
limited resources. Statistical modelling provides one approach to ensuring the effectiveness of inspection and 
maintenance activities. 

Inference for complete systems consisting of thousands of components is a largely overlooked aspect of 
real world inspection planning and inference, due to methodological and computational complexity. Most 
attempts to model a degrading system empirically consider individual system components. Yet characteristics 
of degradation are often common across components, due to component design, age, location, manner of 
operation, etc. A multi-component (or system-wide) model could make use of common between-component 
behaviour to improve the quality of inspection information and hence achieve more efficient inspection 
planning. Inspections at a given time are rarely performed system-wide. A realistic empirical model for a 
large scale system will therefore assume partial inspection data as input. 

Here we develop a model for degradation of a system consisting of multiple dependent components, appli- 
cable to the analysis of irregularly spaced spatially distributed short time series. We make inferences about 
system state per component and model variances, given indirect observations. A Bayes linear approach sim- 
plifies parameter estimation in comparison with a full Bayesian analysis and avoids unrealistic distributional 
assumptions, using model-based simulation. Diagnostic tests assess model fit, and simulation studies based 
on a hypothetical known system are used to evaluate performance. We illustrate the method by modelling 
wall thickness and corrosion rates for corroding pipe-work networks on a full-scale offshore platform, given 
sparse inspection data for component minimum wall thickness. 



Specification of realistic priors and initial values (e.g. O'Hagan et al. ( 2007[ )), for model variances in particular 



(e.g. Farrow ( 2003 ) ) , is problematic in general yet can be highly influential, particularly when data are sparse. 
To avoid over-reliance on poorly-specified priors and initial values, we develop and implement methods of 
inference for the model variance structure. We facilitate variance estimation by making exchangeability 
judgements appropriate for analysis for sparse time-series of irregular partial inspections. 
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Bayes linear methods (Goldstein and WoofF (2007)) applied to dynamic linear models (DLM, Harrison 



and West (19971) offer system- wide modelling of corroding systems using partial inspection data. Little 



et al. (2004a I uses a multivariate DLM to characterise the corrosion of large industrial storage tanks, using 



observations of component minima, and suggests approaches to optimal inspection planning. [Little et al.| 
( 2004b ) describes the application of a spatio-temporal DLM to model the corrosion of an industrial furnace 
using Bayes linear updating. Empirical distance-based estimates for covariances of DLM observation and 
system variances are used, and optimal inspection planning based on heuristic criteria is considered. |Farrow| 
and Goldsteinl (|1993|) discusses Bayes linear methods for grouped multivariate repeated measurement studies 



with application to cross-over trials. Wilkinson (1997) discusses variance learning for a univariate linear 



growth DLM, and Wilkinson and Goldstein (1997) describes Bayes linear covariance matrix adjustment 



for a multivariate constant DLM. Shaw and Goldstein (1999) discusses Bayes linear experimental design 



for grouped multivariate exchangeable systems. Randell et al. (2010) develops a utility-based criterion to 



assess inspection design quality, based on a linear growth DLM with variance learning. Industry guidelines 



(e.g. UK Health and Safety Executive (2002) and ASTM Standard G16) treat the modelling of corrosion 



very generally, yet there is a vast body of engineering literature on this subject. | Zhang and Mahadevan| 
(2000) outlines mathematical expressions for initiation and evolution of different corrosion mechanisms. 



including pitting and cracking. Kallen and van Noortwijk (2005) discusses inspection and maintenance 



decisions based on imperfect inspection within a Bayesian framework, using gamma processes. |Qin and Cui| 
( [2003) considers corrosion of steel structures, and Gasemyr and Natvig ( ,2001 ) presents a Bayesian approach 
using partial inspections only. A number of authors discuss the inclusion of inspection data and expert 



judgement within a risk-based inspection framework. For example, Farber and Sorensen (20001 presents an 



approach to estimating system condition for inspection planning purposes using a combination of inspection 



observations and expert judgement, and Straub (20041 describes generic approaches to risk-based inspection 



of steel structures. Kuniewski et al. (2009) presents a method for compliance sampling. Many approaches to 
corrosion modelling make use of methods associated with extreme value analysis (e.g. | Coles ( 200l| ); Scarf 



(.e.g. 
(e.g. 



and Laycock (1996) provides a review. Glegola (2007) applies the generalised extreme value distribution 



Kotz and Nadarajah (20001), and Lopez de la Cruz et al. (2008) the non-homogeneous Poisson process 
Daley and Vere- Jones |2003[ )). 



A number of aspects of the current work are novel, including Bayes linear variance learning using partial 
observations and Mahalanobis learning for local corrosion variances. The statistical model adopted admits 
a non-linear observation equation. We also believe that the application of the methodology to a full-scale 
industrial system rather than individual components is particularly interesting and informative. 

The article is presented as follows. We start by outlining the motivating application, a full-scale offshore 
platform, in section[2j and Bayes linear methods in section[3j Section|4]introduces the corrosion model and in 
particular discusses the exchangeability judgements made to accommodate irregular partial inspection data. 
Section [5] discusses Bayes linear inference, and the estimation of model variances, involving a Mahalanobis 
distance fitting procedure for local corrosion variance. Estimation and performance of the variance learning, 
in application to simulated data with the same inspection design and covariance structure as the historical 
inspection data, is also discussed. To facilitate future application of the method, section[6]presents a stepwise 
modelling procedure. The model is applied to historical data from the offshore platform in section [7[ and 
diagnostic tools for quality of model fit illustrated. Section |8] discusses findings and proposes possible future 
generalisations. Appendix A| outlines calculations used in computing the adjusted expectations of model 
variances. [Appendix B provides prior values used in the analysis of both historical and simulated inspection 
data. 



2. Motivating application 

We consider inspection of a full-scale offshore platform as a motivating application. For inspection and 
maintenance purposes, the installation is considered as a set of corrosion circuits, each consisting of multi- 
ple components, for inspection. For the current application, we model a system of four corrosion circuits 
consisting of a total of 64 pipe-work weld components. A corrosion circuit typically exhibits long sequences 
of connected components, often with side branches of components. The corrosion behaviour of components 



is influenced in part by their connectivity within a corrosion circuit (see section 4.2). A corroding system 
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may exhibit multiple corrosion circuits with similar corrosion behaviour. This application was discussed 



previously by Randall et al. (20101. 



When components are inspected, inspectors are most concerned with most vulnerable components since 
these are critical to assessing risk of system failure. Measurements are made using non-destructive (e.g. 
ultrasonic) inspection. The use of non-destructive testing means that operations can continue with little or 
no system interruption. The measurement device reports an estimate for the minimum wall thickness over 
a region sometimes referred to as the footprint of the device. Historical data for component minimum wall 
thickness, obtained during inspection campaigns for the period 1998 - 2005, are available. Based on the 
frequency of observations and the requirements for inspection planning, we select a monthly time increment 
for modelling. The sample period therefore consists of 83 time points. 



The actual historical inspection design is given in figure |B.1[ It is clear that inspections are typically 
incomplete and irregularly spaced in time. A total of 174 observations of the system are available. The 
application is used in two ways in what follows. Firstly, it is used to evaluate model performance by 
generating simulated samples using the actual inspection scheme and realistic parameter values. Secondly, 
in section[7j estimates for system parameters based on actual historical inspection data, and system forecasts, 
are produced. For both simulation and historical data analyses, the actual inspection design matrix is used, 
and identical prior specifications made, so that the analyses are comparable. 

Estimates for prior values are obtained from two sources, namely auxiliary historical data from other corrosion 
circuits for the same offshore installation, and judgements made by experienced inspection engineers. (A full 



list of prior values is given in section Appendix B ) 



[Figure 1 about here.] 

3. Theoretical background 

3.1. Bayes linear analysis 

For a complex system, it can be difficult or impractical to make full prior belief specifications. Bayes linear 
analysis allows us to specify and update partial aspects of our beliefs. Bayes linear analysis also provides 
a computationally efficient method for updating beliefs for problems where a full Bayes approach would be 
too difficult or time consuming. It can also be viewed as a generalisation of the full Bayes approach which 
relaxes the requirement for full probabilistic prior specifications. In Bayes linear analysis, expectation rather 
than probability is treated as a primitive quantity; prior beliefs are specified in terms of means, variances and 
covariances. A detailed explanation is given in Goldstein and Wooff ( 2007[). Discussion of the application of 



these methods for the analysis of simulators for large systems is given in Craig et al. ( 2001 ) and Goldstein 



and Rougier (2006) 



Given a vector of data D, the adjusted expectation E£i{B) for a vector B is given by: 

Ed{B) = E{B) + Cov{B, D){V'AT{D))^d - E{D)) 

where the matrix Var(Z3)^ is the inverse of Var(Z?), if invertible, or a generalised inverse otherwise. The 
variance resolved by adjustment, Rvar£)(_B), is given by: 

RvavDiB) = Cov{B, D)(Ya.T{D))''Cov{D, B) (1) 

and the adjusted variance, Varc(_B), by: 

Ywd{B) = Var(B) - Rvarz5(B) 

3.2. Exchangeability and the representation theorem 

The concept of exchangeable events is a crucial component of the subjective theory of probability. In essence, 
exchangeability judgements in a subjective analysis can be used to underpin the types of independence 
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assumptions made in classical inference (de Finetti 1974 1. For Bayes linear analysis, where only partial 



beliefs need to be specified, we can restrict our assumptions for the error structures to exchangeability of 
the first and second order quantities. 

The means, variance and covariances of a second order exchangeable sequence, X = Xi,X2,... , are in- 
variant under permutation. If we assume second-order exchangeability, we can use the the second order 
exchangeability representation theorem, (Goldstein, 1986| ) to express the quantities Xj in the sequence in 



terms of the sum of two random quantities M{X) and TZi{X) which may be viewed as analogous to an 
underlying population mean and discrepancies from the mean respectively, as follows. Given a collection of 
vectors, X = Xi,X2, . . ., an infinitely second order exchangeable sequence with: 

EiX,) = fix, Var(X,) = Ex and Cov(X„ X,) = F^ ^ ^ j 

we can express each Xi as: 

x,^M{x) + n,ix) 

where Ai{X) is a random vector known as the population mean with: 

E{M{X))^fix Var(X(X)) =.Fx (2) 
and the discrepancies TZi{X), themselves second order exchangeable, with: 

E{n,{X)) = and Var(7^,(J^)) = Sx - Fx (3) 
Each pair TZi and TZj is uncorrelated {i ^ j) and each TZi is uncorrelated with Ai{X). 

3.3. Bayes linear inference and its analogy to "full" Bayes 

Bayes linear inference can be understood by analogy to the usual "full" Bayes inference as follows. Suppose 
we observe a quantity D and use it to make inferences about an unknown quantity B expressed as B — A4+TZ, 
where Ai and TZ are unknown random quantities with full prior distribution f{A4,TZ). The prior is related 
to the full posterior distribution f{Ai,TZ\D) by Bayes theorem: 

f{M,n\D) ^ f{D\M,n)f{M,n) 

where f{D\A4,TZ) is the likelihood for D. In the Bayes linear formalism, the second order exchangeability 
assumptions and associated representation theorem B = Ai + TZ provide a structure for B. Bayes linear 
adjustment is then analogous to estimation of the full posterior distribution for Ai and TZ. The specific 
properties of A4 and TZ imposed allow inferences to be made relatively straightforwardly. The non random 
quantities E{Ai) and E{TZ), primitives in Bayes linear inference, are analogous to prior distributions. The 
non random adjusted expectations En^Ai) and E£i{Ai) are analogous to posterior distributions. A fuller 
account is given by jGoldstein and Wooff] ( 2007 section 3.5) 



4. Model 

General framework 

The general framework for our analysis is as follows. We assume that the system to be inspected can be 
partitioned into a set of C components, indexed by c, whose characteristics evolve over T time points, indexed 
by t. We seek inferences about the true system state vector Zid over L locations indexed by / within each 
component. We separate global aspects which affect the whole component from local aspects. This allows 
us to distinguish between different model characteristics. 

Global effects 

The global effects model captures the most important features of, and relationships between components. 
Global effects evolve in time as a dynamic linear model (DLM) with system evolution matrix G: 

Qt = GQt-i + eet (4) 
where Qt and egt are vectors over components with elements Oct and tgct respectively. 
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Local effects 

The local effects model describes spatial variability in detail. Local effects rid evolve in time for some 
function g as: 

rict = g{ric(t-i)) + f^rict (5) 



True System State 

We model Zid as the sum of global and local effects: 

Zlct = FcOt + rict 

where Fc is the cth row of the matrix F of linear combinations of global effects parameter vector Ot- 



(6) 



Observations 

For each component and time, we can choose at a cost to observe a function / of the true system state vector 
over locations Z^t — ( ■ ■ ■ Zict • ■ • ) with error: 



^ct — f {Zct + f'Yct) 



(7) 



where eyct is a vector of measurement errors across locations. We assume that the function / is non-linear 
and separable in the following sense. For vectors a and 1 (a vectors ones) of the same length, and a scalar 
&, / can be decomposed as: 

/(a + 61) -/(a) +6 

Thus: 

/ {Zct + f-Yct) = / (^"6* + rict + evct) 
= FQt + f {rict + evct) 

We note that separability applies to many functions for summarising data, including the mean, the median, 
maximum, minimum and other quantiles. 

Complete Model 



Observation Equation 
True System State 
Global Effects Model 
Local Effects Model 



Yet ^ f {Zct+ ^Yct) 

Zlct = FcQt + rict 

Qt = GQt-i + est 

rict = g{ric{t-i)) + erict 



We assume that errors egt, ^rict and eyict are mutually uncorrelated in time. The elements egct of egt are 
correlated across components with Var(eect) = Sg, but the trict are independent across components. We 
assume Yai{erict) = and Yai{tYict) = both scalar constants. Specification of error structure in terms 



of second order exchangeability is discussed in section 4.2 



4-.2. Exchangeability judgements 

Non-linearity of the observation equation makes formal Bayesian inference difficult. Coupled with the large 
system size and resulting difficulty of prior specification, full Bayesian calculations become intractable. This 
is particularly important when considering inspection design requiring fast evaluation since many samples 
are needed to assess each of a large number of design choices. Here, we proceed using Bayes linear methods 



(section 3.1 ) in conjunction with plausible exchangeability judgements. The latter provide a means to define 



model variances from squared linear combinations of observations since they provide access to model evolution 



errors (section 5.2). In order to invoke second order exchangeability we view the observations in time and 
space as part an infinite exchangeable sequence. We make the following exchangeability judgements: 
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Judgement: Second order exchangeability of squared errors e^^^ over time for each component. 
This leads to representation statements for the squared residuals of every component: 

4ct = yect^M{Vg,)+'Rt{Vec) 
which can then be further decomposed. 

Judgement: Second order exchangeability of AiiVgc) across components such that: 

M{Vgc) = Wee = M{We)+n,{We) 

where: 

E{Wec) = ti-We , Var(M^ec) = ^Wa and Cov(PFec, Wed) = ^We for c^c' 

so that: 

E{el^^) = E{Vect) = E{M{Vec)) = E{Wec) = t^w, and Var(A^(F,,)) = Eh-, 
Then, from equations [2] and [3] respectively: 

YaiiMiWe)) = r^, and Var(7^c(VKe)) = T.Wo - ^We 

and for each e^^j we can write: 

el,=M{Vec)+nt{Vg,) 

= M{We)+nc{Wg)+ntiVec) (8) 

Judgement: Covariance between residuals 

We might try to learn about the global effects covariance matrix Sg in full generality. However in practice 
this proves difficult. Instead we choose to express Sg in terms of known correlation Hg and unknown second 
order exchangeable variances A4{Wg) and TZdWg). The elements Sg^^, of the global effects covariance matrix 
Sg are given by: 

Se^.> = Wglwlllg^d 

= iMiWg)+nciWg))'HMiWe)+nc'iWg))mg,d (9) 

where 

/ N / Ilgcc' for t = t' 

corr(6,,„e,,,,0 = | for 

The prior expectation E(Sg ,) — T.g , . To estimate Sg using a sample, we must estimate A4{Wg) and 
Ti-ciWg). In practice, where it is too difficult to learn about TZdWg) we also treat it as known. To update 
beliefs about Sg, we therefore only need to learn about Ai{Wg). 



4-.3. Example: Corrosion framework 



For the motivating corrosion application (section [2]), ultrasonic and radiographic inspection typically gen- 
erates values for the minimum wall thickness (or maximum pit depth) corresponding to the area inspected 
(the "inspection footprint" ) . We therefore adopt the specific model framework: 



Observation Process: 

True System State: 
General Corrosion Model: 

Local Corrosion Model: 



Yet = mm(Zct + eYct) 

Zict — Xct + rict 
Xt = Xt-i + a* + ext 
at = "(i-i) + Eq* 

"Tlct = nc(t-l) + (-rlct 



(10) 

(11) 

(12) 
(13) 
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where Zd — (Zict) and eyct are vectors over locations. With reference to section 4.1 / is the minimum 



function over set L, g is the identity function, and functions F, G and parameter vector 9t take the form: 



_ <^t J \ Oc Ic 

where Ic and Oc are C xC identity and zero matrices respectively, and is a vector of wall thicknesses over 
components. We seek inferences about true system state Zi^t and corrosion rate vector at over components. 
The error terms exct, ^act and e^ict control evolution of system level, slope and local level adjustment 
respectively. 



We make the analogous exchangeability assumptions as in section [4~2| so that: 

M{Wx)+ Tl,{Wx) + TltiVxc) 

as given in equation [s] In principle it is possible to learn about e^^j and e^^^ separately. In practice, given 
sparsity of data, this can be difficult. For simplicity we fix the ratio of mean variances: 

M{Wa)) = \M{Wx) (14) 

The correlation structure of wall thickness evolution error IIxcc' is assumed to take the form of a linear 
combination of three terms, reflecting an underlying universal correlation po between all pairs of components 
(regardless of the circuit (s) to which they correspond), a circuit correlation pc between all pairs of compo- 
nents within the same circuit, and a correlation p£, which decays exponentially at rate v > Q with distance 
s (measured in terms of the number of intervening components along the circuit between the components) . 
So the covariance between components c and c' is: 

IIxcc' = Po + Pc^cc + PDe-^""""' (15) 

where 5cc' = 1 if components c and c' are in the same circuit, and Scc' is the distance between the components. 
For simplicity we assume 11^ ~ lix and then: 





n„ 



such that Sx and are defined using equation [9] The correlation matrix IIx for the offshore application 
is illustrated in figure [B?2} 

[Figure 2 about here.] 



5. Inference 

The objective of the analysis is to learn about the characteristics of model parameters in equations 9-12. 
We achieve this using an iterative scheme explained in sections 5.1-5.4. We use Bayes linear mean updating 
in section [5.1 1 to learn about the true system state Z^t and global effects parameter vector 9t assuming that 
all variance parameters in the model are known. In section [5. 2 1 we show how Bayes linear variance learning 
is use to adjust beliefs about global effects model variance parameter M{Wg) assuming known local effects 



variance parameters. In section 5.3 we use Mahalanobis variance learning to estimate local effects variance 
parameters given updated global effects model variance M{Wg). Model fit and prior assumptions are 
assessed using discrepancy diagnostics discussed in section 5.4. 

5.1. Bayes linear mean updating 

Given a collection of observations of one or more system components at one or more times expressed in 
vector form as F, we adjust beliefs about true system state Z""" and global effects level 9t by calculating 
EviZ^n and Eviect): 

EviZ^"') = S(Z™") + Cov(Z™", y)[Var(y)]t(y - E{Y)) (16) 
Eviect) = EiOct) + Cov(0,t,r)[Var(r)]t(y - E{Y)) (17) 
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where = mmi{Zct) over locations. In a similar fashion we can calculate adjusted variances: 



ct ) 



Vary(Z™"') = Var(Z^7") - Cov(Z^f y)[Var(y)]tCov(r, Z, 
Vary(0et) = Var(0,t)-Cov(^^,t,y)[Var(r)]tCov(y,0,t) 

Estimates for (co-)variances Var(y), Var(Zict), Var(0ct), Cov{Zict,Y) and Cov{9ct,Y) for corresponding 
components and times are obtained by simulation under the model using prior beliefs. Details of the sim- 
ulation procedure are given in section |6] Briefly we simulate realisations for the complete system from the 
model given in section [4] given prior specification. Realisations are then use to calculate empirical estimates 
for expectations, variances and covariances used in Bayes linear adjustment. We can use simulation to run 
the model forward in time beyond the period of the data or to predict characteristics of unseen components 
in a straight forward manner. 

5.2. Bayes linear variance learning 

As illustrated in section 4.2, the Bayes linear approach may also be used for variance learning. Here, develop- 
ment of Bayes linear variance learning is considered for the corrosion modelling application for concreteness, 
since the corrosion application illustrates all key modelling features. Other applications will require modi- 
fication of calculation details in general, but the methodology remains applicable. We seek expressions for 
squared residuals e^^^ from sample data corresponding to partial system inspections with which to learn 
about population mean variances M{Wx), and thereby M{Wa) and also. Since the general corrosion 
DLM part of the model is invertible we can take linear combinations of observations to isolate expressions 
for exct 6ven when observations of the system are irregularly spaced in time. Using expressions for e^^j 
thus obtained, we adjust our beliefs about M{Wx) using observed data and assumed known values for 
as explained below. In section [5. 3| we present a fitting procedure using a Mahalanobis distance criterion to 
select an optimal combination of error variances. 

Motivation 

We motivate the theorem below by considering the case of full regularly spaced inspections, with the obser- 
vation equation expressed as: 

Y,t = X,t + M,t (18) 

where: 

Met = min (r/ct + evict) 

consider taking differences of observations in time. Let y^j*-* be the ith step difference for component c at 
time t defined by: 

^ct ^ =Yct — Ycit-i) 



with the analogous definition for M^l\ Then using equations 9 and 18 one and two step differences are seen 
to be: 

ri/) ^ - Fc(t-i) - X,t - -I- Met - M,(t„i) 

= Xc(t-r) + o-ct + exct - Xc(t-i) + M^l^ 
= act + exct + Mrf ^ 



and: 



ac(t-i) + e«ct + £xct + Mj:l^ 



— ^ct — Yc(t~2) — Xct — Xe(t-2) + M^^^ 

(2) 

= '2ac{t-i) + ^act + exc(t-i) + exct + M^^ 
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The linear combination Y^'^'^ — 2Y^l^ isolates the residual error structure except for the term m'^'' — 2M^l^: 

The expectation of the square of the former linear combination is found using the facts that residuals 
are mutually uncorrelated in time, together with the exchangeability assumptions from equatioijS] and the 
assumption of a fixed ratio of population variances in equation |14| Then: 

E{y}^^ - 2Y}pf =E (e^,t + exc(t-i) - exct + M^f^ - 2^1^^^" ' 



--E(e 



--E\y\{M{Wx) +nc{Wx) + UtiVxc)) 

+ M{Wx) + n,{Wx) + nt-i{Vxc) 

+ M{Wx) + n,{Wx) + TltiVxc) + - 2A4l^f 



--E ((A + 2)M{Wx)) + E ((M(f) - 2M^J^r ^ 
--{\ + 2)^w. + E ((Mi') - 2MiP f) (19) 



Using the corresponding result for irregular inspection, outlined in appendix [Appendix A.l we arrive at 
the following theorem. 

Theorem 5.1. Consider incomplete inspection of a system yielding observations {Yct^,Yct^, . . . ,Yct^ } at 
times {ti, i§, . . . ,t^ } for components c = 1, 2, . . . , C. The adjusted expectation E^{M{Wx)) for A4{Wx) 
given A is: 

E^iMiWx)) = E{M{Wx)) + Cov(7W(W^x), A)(Var(A))t(A - E{A)) 
where A is a C x 1 vector over components with elements Ac.' 



i=3 



Further, E{MiWx)) = ^x, Cov(X(M/x), A) = l^iT^ - 2)Twx and: 
p.x ^ ktmkf-l'^)(2Xm^-2Xl'^-6-X) 



+ {I'rfEm!:^})') + {KfEm^lf) - 2kmE{M^}M':^}) 
where k'i — t^ — t^-i, li ~ ^1 ~ ti-2 and Y^^J — Ycf; ~ Yd". with a similar definition for M^'j.c . 



Proof. E{A4{Wx)) = IJ-Wx by definition in section 



4.2 



Derivations of expressions for E{A) and Cov(A^(Wx), A) 



are given in appendix [Appendix A.l] and [Xppendix A.2| □ 

Corollary 5.2. In the case of full regular inspection yielding observations Y^t at times t — 1,2, ... ,r for 
components c, the expression for Ac is: 

t=3 



Proof. Let k'( — 1 and 1^ — 2 for each t"^ in theorem 5.1 and E{A) is found as in equation 19 



□ 
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5.2.1. Corrosion Example 

For known local corrosion variance, E^, we use simulation with prior beliefs to estimate E[{m';i^})% E[{m';^})% 
and E[M^}m^}] and Var(A))^. Details of simulation procedure are given in sectionjoj As illustration, con- 
sider estimating ^A{Wx) using E^{M.{Wx)) for a simulated case using the actual mspection design from 
section 2, with M.{Wx) set to 0.1'^. For each of 50 independent realisations of the systems over all compo- 
nents and times, we calculate A at the inspection times, which is then used to calculate E^{M.{Wx))- The 
mean value of E^{M.{Wx)) is 0.0994^, with empirical 5% and 95% values (from simulation) of respectively 
0.0872-^ and 0.1128^, consistent with the known value. 

This procedure allows variance learning for the general corrosion DLM when local corrosion variances are 
known. In practice this is not the case. Selection of optimal combination of local and general variance 
estimates consistent with data is discussed next. 



5.3. Mahalanobis variance learning 



As for general corrosion in section [5T2| we would ideally use a Bayes linear scheme to update local corrosion 
variances also. However the non-linear nature of the observation equation (equation 9) renders direct estima- 



tion impossible. Instead we adopt a fitting procedure based on Mahalanobis distance (Mahalanobis (1936)), 
exploiting the estimated covariance structure, to estimate combinations of local and general corrosion error 
variances consistent with observational data, given an estimate of measurement error variance. 

For each of a set of p candidate values for local corrosion variance S^, we calculate E^{Ai{Wx)) by Bayes 
linear adjustment. Adopting E^{M{Wx)) as an updated estimate for M{Wx), and XE^{M{Wx)) for 
A4{Wa), we re-simulate to estimate E{Y), Var(y) and the ratio of Mahalanobis distance to its expected 
values, termed the discrepancy ratio, H: 

^ ^ [Y - E{Y)YY'^v{Y)\Y - E{Y)) 

rank(Var(r)) ^ ' 

as outlined in section [6j We select the candidate which yields a discrepancy ratio nearest to its expected 
value of unity. 

[Figure 3 about here.] 
[Figure 4 about here.] 

We illustrate joint Bayes linear and Mahalanobis variance learning for the simulated example from section 
5.2.1 in figures [K3| and rB.4[ Figure [b73| shows Bayes linear variance learning for wall thickness variance 
E^{AA{Wx)) as a function of local corrosion variance E^- The simulated example has true values oi M.{Wx) 



and both 0.1^. Figure B.4 shows the discrepancy ratio H as a. function of local corrosion variance 
for the same example. Comparison with the expected value of 77, shown as a horizontal line, suggests the 
estimate E^ = 0.08^. Figure B.3 then provides the corresponding estimate E^{J^{Wx)) = 0.11^. 



5.4. Diagnostics 

Diagnostics are essential to confirm the adequacy of model fit. Using Mahalanobis distance, an adjust- 



ment discrepancy Disy(X) (see page 122 of [Goldstein and Wooff (2007)) for adjusted expectations can be 
calculated: 

V)isy{X) = {Ey{X) - E{X)yRYa.TYX-\EYiX) - E{X)) 

where X is a vector of wall thicknesses for all components and times and RVaryX is defined in equation [T] 
Analogously, Disr(-^) and Disy(a) can be calculated for vectors Z and a over all components and times. 



We can also compute discrepancy directly on the data (see page 100 of Goldstein and Wooff (2007)): 



Dis(r) = (r - E{Y)fYaT{Y)-\Y - E{Y)) (21) 
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The expected value of both Disy (X) and Dis(F) is unity. Forms for both discrepancies can also be evaluated 
for subsets of variables, useful, e.g., to spot individual outliers. Discrepancies are calculated on the prior 
specification, adjusted expectations of global effects terms and, subsequent to learning about local effect 
variance parameters, quantifying model adequacy at each stage. General thresholds for discrepancy measures 



do not exist although as a heuristic the three sigma rule (Pukelsheim 1994 1 states that for any uni-niodal 
continuous random quantity, X, P{\X - E{X)\ < 3^Var(X)) > 0.95. 



[Figure 5 about here.] 



Figure B.5 shows Dis(y) for a realisation of the simulated data discussed in sections 5.2 and 5.3. The 
expected value of Dis(y) is unity, shown as a horizontal line. In this realisation, only one observation exceeds 
the warning limit suggesting that prior specification and data are consistent. More general diagnostics for 
variance learning and overall model fit can be specified similarly. 



6. Inference procedure 



The procedure to apply the corrosion model (section 4.3 ) is explained below. Firstly we make prior specifi- 
cation and carry out simple diagnostic checks to confirm consistency of the data and priors. We then update 
model (co-)variances, before updating means. Finally we calculate diagnostics to assess model fit. 

To simulate from the model, distributional assumptions are necessary below (indicated using notation iid(.). 
It is advised that a variety of different distributional forms be examined to enable an informed choice 
about (co-)variance specification. Plausible distribution will generally lead to robust estimates of low order 
moments, even if full probability specification would not be appropriate using these distributions. 

1. Specify prior and starting values to be examined 

• wall thickness variances prior, ^Wxj ^vFx ^-^d (see section 4.2), and Hx ^, (see equation 15) 

• corrosion variance variance prior. A, Hq,^^, (see section 4.3 ) 

• measurement error variance Ey (see section 4.1) 



• candidate values for p local corrosion variances to consider (see section 4.1) 

• starting values over components for Xt and a* at t = (see section 4.3 ) 



2. For each of p choices of local corrosion variance (see section 4.1 ) 



(a) Evaluate variance matrices using simulation under prior specification 

• simulate Sx from ud(^vi'x i i i n^^^, ) 

• simulate Sa from iid(A/ivi/x ; ^^Wx ) ^^Wx i ^a^^, ) 

(b) Run model specified in section 4.3 forward from t = 

• ea.t ^ iid(0,5a) 

• a.t = a.t-i + Cq.* 

• ex.t ~ iid(0,5x) 

• X,t = X,t_i + a,t + ex.t 

• er..t ~ iid(0, T,r) 

• r,,t = r..t_i + er..t 

• ey..t iid(0, Sy) 

• y.t = min/ {X,t + r..t - 1 + ey..*) 

(c) From simulation calculate sample estimates E{Y) and Var(y) over N realisations and check 
consistency with prior by computing discrepancy Dis(y) (see equation [2l|) 

(d) Calculate A and Var(A) over N realisations and calculate E^{M.{Wx)) for each (see theorem 



5.1) 



3. Select optimal from the set of p candidates as that which yields the discrepancy ratio (see equation 



201 nearest to unity 



4. Simulate to update wall thicknesses Xt and corrosion rates at for any t of interest 
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repeat step 2a and 2b using optimal and fiwx = E^{M{Wx)) for N realisations 

calculate sample estimates from simulation for £'(Z"""), E{act), Cov(Z™"\ Y),Cov{act,Y), Var(y) 
and E{Y) 

calculate adjusted expectations from true wall thickness and corrosion EY{oict) (see 

section 5.1 1, for components and times of interest. This could include forecasting forward in time 



or prediction of components not observed 



5. Calculate discrepancy ratio H (see equation 20) for estimated model variables to confirm consistency. 



7. Application 



Having illustrated model performance for simulated data, we now apply it to actual historical inspection 
data for the offshore platform, using the procedure described in section [6} T he first step of the analysis is 
specification of prior values. Recalling that the inspection design (figure B.l I and system prior specification 
(section [2) corresponding to the historical data were also used for simulated data, results equivalent to those 
in figures B.3 B.4 and B.5 are shown in figures B.6 B.7 and B.< 



below. 



We assess the prior discrepancy of observations as described in section |5.4[ To achieve this, realisations 
under the model are generated as described in steps [2a| and |2b| of section [6] Discrepancies are shown in figure 
B.6 to be compared with figure B.5 for a single system realisation. With one exception, historical inspection 
data are consistent with prior expectations. At time point 11, corresponding to the first observations of 
components 35 — 51, all in a particular corrosion circuit, high values of discrepancy may refiect poor prior 
specification for that circuit at that time. 



[Figure 6 about here.] 



Starting values for wall thickness Xq and corrosion ao were specified using historical data. We evaluate 
adjusted expectation for wall thickness variance E^{Ai{Wx)) as described in theorem |5.1| for different 
local corrosion variance using variance learning (section 5.2). Results are given in figures B.7 and B.8 
(analogous to figures B.3 and B.4 for the simulated data). The value of 0.1^ for E^{Ai{Wx)), to be used as 



an estimate for prior wall thickness variance M{Wx) is the same as its prior value (Appendix B ). However, 
the updated value of 0.4^ for is considerably larger than its prior (0.1^). Wall thickness forecasts and 
inspection designs derived using updated estimates for corrosion variances will be different in general to 
those based on constant (prior) variances. 



[Figure 7 about here.] 



[Figure 8 about here.] 



[Figure 9 about here.] 



The forward model together with parameter estimates from historical inspection data allows forecasting of 
future system state. Figure [B79| illustrates the effect of incorporating variance learning for a single component 
of the offshore platform. The critical wall thickness corresponding to component failure is shown as a hori- 
zontal dotted line at 4mm. Actual inspections of the component are shown as black circles. The grey shaded 
area corresponds to a 95% credible interval for underlying wall thickness based on simulation using prior 
beliefs (with mean shown as a white line). Light grey broken lines correspond to the adjusted expectation 
i?y(^™™) and 95% credible interval for underlying wall thickness based on Bayes linear adjusted variance 
Vary(Z^™). Solid black lines correspond to the same adjusted values following a variance learning step. 
Since variance learning inflates local corrosion variance, credible intervals incorporating variance learning are 
wider. The rate of wall thickness loss, and the expected time to crossing of critical wall thickness are also 
affected by the incorporation of variance learning. Notwithstanding large uncertainties, in this application, 
variance learning suggests that component life may be approximately 20 months shorter than otherwise 
anticipated. 
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Figure B.IO illustrates the effect of incorporating variance learning for 3 components of the offshore platform, 
one of which is not directly observed. The first component is the same component as shown in figure |B.9[ 
The grey shaded area corresponds to a 95% uncertainty interval for underlying wall thickness based on 
simulation using prior beliefs (with mean shown as a white line); these are the same for each component. 
Light grey broken lines correspond to the adjusted expectations £'y(Z™'") with 95% uncertainty intervals 
based on Bayes linear adjusted variance VarY(Z™™). Solid black lines correspond to the same adjusted 
values following a variance learning step. The effect of the correlation between components can be seen 
in changes exhibited by adjusted wall thickness for components at times with no inspections. The third 
component shows the effect on a component which is not directly observed. In the case of mean updating 
(without variance learning), the corrosion rate is reduced, increasing predicted remnant life. In the case of 
variance updating and mean updating, the local corrosion rate is increased and the corrosion rate is reduced. 
The net effect is to decrease predicted remnant life. 

[Figure 10 about here.] 



8. Discussion 



We present a general framework for modelling the evolution of a large scale physical system consisting of 
multiple dependent components. Sample data consist of short time series of irregularly-timed partial in- 
spections across components. Realistic modelling requires good specification of priors and initial values, 
particularly if data are sparse. Simulation and Bayes linear analysis provide knowledge of the full system 
given partial prior specification, considerably more simply than full Bayesian inference. Given second order 
exchangeability judgements, variance structure learning for short time series is achieved using linear combi- 
nations of observations to access squared residuals where possible. When the observation equation involves 
a non-linear function a Mahalanobis fitting procedure can be used estimate local effects variance. 

For the model form given in section [4j once distributional assumptions are made, it is straight forward to 
simulate realisations to estimate means, variances and covariances. 

Second order exchangeability judgements used are weaker than independence assumptions typically made in 
classical inference. Exchangeability across components allows us to exploit large numbers of components for 
variance structure learning, particularly advantageous when data are sparse. Diagnostic tests assess model 
fit and validity of exchangeability judgements. Specific knowledge of system characteristics allow a more 
detailed prior specification, e.g. partitioning the system into sets of exchangeable components, or specifying 
a parametric form for certain variances in time. With sufficient data, between-component exchangeability 
judgements may not be necessary. Separate inferences about variances of individual components could be 
made. We illustrate the method in application to corrosion modelling of a full scale offshore facility. Forecasts 
of remnant life and its uncertainty incorporating variance learning are shown to be materially different to 
those obtained without variance learning. 

Various possibilities for generalisations of the model exist. For example, the corrosion model ignores cor- 
rosion initiation, which could be incorporated using an appropriate (e.g point process) representation. For 
applications other than to corrosion, alternatives to the (separable) minimum function in the observation 
equation and DLM structure will require tailored linear combinations to access squared residuals required 
for variance structure learning. Nevertheless, the strategy of Bayes linear variance learning, where possible, 
supplemented by Mahalanobis learning, will be applicable. 
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Appendix A. Calculations for Er){M.{Wx)) 



Appendix A.l. £'[A] 

The general corrosion DLM (equation |4| can be rewritten: 



X 

a 



1 k 
1 



X 

a 



fc-i 

c(t-fc) i=o 



1 i 
1 



then 



(10) X,t. -(10) ^c(t^-fc,^) + M^} 

( fcf ) Xc(tJ-iJ) + ( /Cf ) e,(t._fcc,tc_,c) 

(1 ) 



where: 



then similarly: 



fc?-i 



Cc(t=,t=-fcf) — E 



1 i 

1 



=(tf-i) 



yif - ( 1 ) x,,c -(10) + M^fi 

= ( ) Xc(4c_;c) + ( 1 fcf ) ^c(tf-fe,^t=-/J) 

+ (10) Cc(t^t^-;cJ) + M(f/ 



(2) 

Now to eliminate the effects of the wall thickness term we do fcf F^^c — 



(1). 



jLC-^(2) _ 7C-i^(l) _ 7,C 



( Zj'' ) Xc(te_;c) + ( 1 fcf ) S^cdl-k-r.tl-lD 

10)^, 



c(t-i) 



( fcf ) + ( fcf ) C,(tc_,c_,c_;c) 

+ (10 )^,(,c,,c_,.)+AfWl 



We have that: 

" I ij^ j 
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So: 



Then: 



2=0 

E 



1 i 

1 



2=0 

E 



1 i 
1 



Var 



1 i 
1 



ex + Cc 







( 








1(1 ; 







kf + hkfikf + l){2kf + 1)A ifcf (fcf + 1)A 



fcf (fcf + 1)A 



kfX 



(l)^2l 



+2(fcr) 



h ~ ki + 



6 

AQ? - fcf)a? - fcf + 1)(2G? - fc^ + 1) 
6 

\in - kf){n ^ k-r + 1) 



+ {kif{ki - l^fHi - ki)X^Wx 
+ {l^^fE[{Mil}S] + m'E[{M^fJf 
_ k^ViiM - lf)C2X{ktf - 2X11 - 6 - A) 



2k1llE[M^jM^^l] 



flWx 



Appendix A. 2. Cov[A^(Wx), A] 

CoY[MiWx),A] = {Cov[M{Wx),Ai],Cov[M{Wx),A2],...,Cov[M{Wx),A^]) 



Cov[M{Wx),A,] - E i^Cov [M{Wx), {k^Y^ - VrY^ 



i=3 



from [Appendix A.l| we know that: 



(l)^2 



( {ki-ll) ) ^c(t<^,t^-fc^'^) + ( ki ki{ki-l'i) ) ^c{tf-fe<=,t^-if) 



+ klM'^ - l^M, 



(1) 



which means: 



Gov 



[m(Wx), {KY^ti - ny^l^f] = Gov \m(Wx), (( (kl - 11) ) C,(tc_,o_fee))' 



■ Gov 
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we will consider: 
A: 

B: 

A: 



Gov 



M{Wx),{{ (kf-m )^,(,c,,c_fec))^ 



Gov 



M{Wx),{{ kt KKM-n) )L(ti-kt^-it)y 



B: 



Gov 



M{Wx),{{ (K-ID )?e(t=,t=-/c=))' 



: Gov 



J=0 



: Gov 



1 



<^X + So 



j=0 

■-{kl-llf E Cov[>l(Wx),>l(Vya)+7?.c(H'a)+7et(yac)] 

-(K-^f E (i + l)'Cov[A1(iy^),A^(iyc«))] 



^OiCtf—j 



So: 



Gov 



: Gov 



: Gov 



\ 3=0 ^ 



^X + 



ct9 — k^ —j 



M(Wx), r • • / 



- ifcf )(fe?)2r,y^ + Gov 



.=0 - I1)el^^^_ 



{It - ktmyvw^ + 



HKri^ - Km - K + mm - K) + i)Twx 



^mf(ii - Km - K + 1)^^? - n)^wx 



+x{kt?{kt-nf{it-kt)Twx 



Gov 



M{Wx),{k'rY^:J-l'iY^J) 



= A. + B. 

klim - l^mkt? - - 6 - X)Tw. 
6 
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Therefore: 



Cov[M{Wx),A,] = E^r^^x 

1=3 * 

= {Tc — '2)Twx 

and: 

Cov[M{Wx),A] = ( 

^Wx ^Wx ■ ■ ■ ^Wx ) 

Appendix B. Prior values for offshore structure application 



number of components 


N 


64 


number of time points 


T 


83 


total number of inspections 




174 


minimum level of correlation between components 


Po 


0.2 


distance effect scaling parameter 


PD 


0.3 


circuit correlation scaling parameter 


PC 


0.5 


ratio local to general corrosion 


A 


0.02 


measurement error variance 




0.16^ 


prior mean for M.{Wx) 


P-Wx 


0.1^ 


prior variance for A4(Wx) 


^Wx 


le-3 


prior covariance for M.{Wx) 


Fwx 


5e-4 
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List of Figures 



B.l Inspection design for the offshore apphcation, consisting of 64 components over 83 time points. 
Black hnes correspond to 174 observations of the system 



B.2 Correlation matrix Hx for the offshore application. The 4 blocks present correspond to the 
fom' corrosion circuits 

B.3 Bayes linear variance learning for wall thickness variance E^{A4{Wx)) as a function of local 
corrosion variance, S,., for simulated data, shown on square- root scale. The true values of 
■M{Wx) (for which E^{Ai{Wx)) is an estimate) and are both 0.1^, as shown by the 
dashed horizontal and vertical lines. The mean estimate for E^{A4{Wx)y^^ is shown as a 
solid line, and the shaded region corresponds to a 90% uncertainty band for _E^(A^(Wx))^^^ 
bounded by the 5th and 95th percentiles derived from simulation 



B.4 Discrepancy ratio, H, as a function of local corrosion variance for simulated data. The 
mean estimate for H is shown as a solid line, and the shaded region corresponds to a 90% 
uncertainty band for H bounded by the 5th and 95th percentiles derived from simulation. The 
expected value of H shown as a horizontal line is 1, suggesting = 0.08^ and corresponding 



E^{A4{Wx)) ~ 0.11'^, (from figure B.3 1. Vertical lines indicate the mean (solid) and 5th and 
95th percentiles (dashed) for the particular choice of Tir in individual realisations. The true 
values of M{Wx) and are both 0.1^ 



B.5 Component-wise discrepancy Dis(y) for a typical realisation from the simulated data. The 
expected value of Dis(y) is unity, shown as a horizontal line. Also shown is the horizontal 
line corresponding to |1 — Dis(K)| = 3, serving as a warning limit for unusually large values 
of discrepancy. In this realisation, only one observation exceeds the warning limit 



B.6 Component-wise discrepancy Dis(y) for the historical inspection data. The expected value of 
Dis(y) is unity, shown as a horizontal line. Also shown is the horizontal line corresponding to 
|1 — Dis(y)| = 3, serving as a warning limit for unusually large values of discrepancy. At time 
point 11, corresponding to the first observations of components 35 — 51, all in a particular 
corrosion circuit, high values of discrepancy may reflect poor prior specification for that circuit 
at that time point 

B.7 Bayes linear variance learning. E^{A4{Wx)) as a function of local variance for the histor- 
ical inspection data, shown on square-root scale 



B.8 Discrepancy ratio, H , as a, function of local corrosion variance for the historical inspection 
data. Figure suggests local corrosion variance S^. = 0.4^ and corresponding adjusted general 



corrosion variance E^{J\A{Wx)) = 0.1 (from figure B.7) 



B.9 Effect of incorporating variance learning for a single component of the offshore platform. 
Critical wall thickness is shown as a horizontal dotted line at 4mm. Actual inspections of the 
component are shown as black circles. Grey shaded area corresponds to a 95% uncertainty 
bands from simulation for underlying wall thickness based on prior beliefs (with mean shown 
as a white line). Light grey broken lines correspond to adjusted expectation i?y(Z^™) and 
95% uncertainty bands for underlying wall thickness based on Bayes linear adjusted variance 
Vary(Z™™) without variance learning. Solid black lines correspond to the same adjusted 
values following a variance learning step 



B.IO Effect of incorporating variance learning for 3 components of the offshore platform, one of 
which is not directly observed. Critical wall thickness is shown as a horizontal dotted line 
at 4mm. Actual inspections of components are shown as black stars. Grey shaded areas 
corresponds to 95% uncertainty bands from simulation for underlying wall thickness based 
on prior beliefs (with mean shown as a white line). Light grey broken lines correspond to 
adjusted expectations i?y(Z^") with 95% uncertainty bands based on Bayes linear adjusted 
variance Vary(Z™*") without variance learning. Solid black lines correspond to the same 
adjusted values following a variance learning step 
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Figure B.l: Inspection design for the offshore application, consisting of 64 components over 83 time points, 
correspond to 174 observations of the system. 
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Figure B.2: Correlation matrix IIx for the offshore application. The 4 blocks present correspond to the four corrosion circuits. 
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Figure B.3: Bayes linear variance learning for wall thickness variance E^{Ai{Wx)) as a function of local corrosion variance, 
Sr, for simulated data, shown on square-root scale. The true values of Ai{Wx) (for which E^{M{Wx)) is an estimate) and 
Sr are both 0.1^, as shown by the dashed horizontal and vertical lines. The mean estimate for E^(Ai{Wx))^^^ is shown as 
a solid line, and the shaded region corresponds to a 90% uncertainty band for E^(M[Wx))^^'^ bounded by the 5th and 95th 
percentiles derived from simulation. 
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Figure B.4: Discrepancy ratio, H, as a function of local corrosion variance for simulated data. The mean estimate for H 
is shown as a solid line, and the shaded region corresponds to a 90% uncertainty band for H bounded by the 5th and 95th 
percentiles derived from simulation. The expecte d va lue of H shown as a horizontal line is 1, suggesting Sr = 0.08^ and 
corresponding _E^(yV((Wx)) = 0.11'^, (from figure B.3l. Vertical lines indicate the mean (solid) and 5th and 95th percentiles 
(dashed) for the particular choice of Sr in individual realisations. The true values oi M(Wx) and Sr are both 0.1^. 
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Figure B.5: Component-wise discrepancy Dis(y ) for a typical realisation from the simulated data. The expected value of Dis(y) 
is unity, shown as a horizontal line. Also shown is the horizontal line corresponding to |1 — Dis(y)| = 3, serving as a warning 
limit for unusually large values of discrepancy. In this realisation, only one observation exceeds the warning limit. 
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Figure B.6: Component-wise discrepancy Dis(Y) for the historical inspection data. The expected value of Dis(y) is unity, 
shown as a horizontal line. Also shown is the horizontal line corresponding to |1 — Dis(y)| = 3, serving as a warning limit for 

unusually large values of disercpaney. At time point 11, corresponding to the first observations of components 35 — 51, all in a 
particular corrosion circuit, high values of discrepancy may reflect poor prior speciflcation for that circuit at that time point. 
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Figure B.9: Effect of incorporating variance learning for a single component of the offshore platform. Critical wall thickness 
is shown as a horizontal dotted line at 4mm. Actual inspections of the component are shown as black circles. Grey shaded 
area corresponds to a 95% uncertainty bands from simulation for underlying wall thickness based on prior beliefs (with mean 
shown as a white line). Light grey broken lines correspond to adjusted expectation and 95% uncertainty bands for 

underlying wall thickness based on Bayes lineaj: adjusted variance Vary(Z5'") without variance learning. Solid black lines 
correspond to the same adjusted values following a variance learning step. 
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Figure B.IO: Effect of incorporating variance learning for 3 components of tiic offshore platform, one of which is not directly 
observed. Critical wall thickness is shown as a horizontal dotted line at 4mm. Actual inspections of components are shown 
as black stars. Grey shaded areas corresponds to 95% uncertainty bands from simulation for underlying wall thickness based 
on prior beliefs (with mean shown as a white line). Light grey broken lines correspond to adjusted expectations Ey(Z^'^) 
with 95% uncertainty bands based on Bayes linear adjusted variance Vary(ZJJ'") without vajriance learning. Solid black lines 
correspond to the same adjusted values following a variance learning step. 
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